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Abstract 

O 

■ We study the quantum- mechanical evolution of the nonrelativistic oscillator, 

. rapidly moving in the media with the random vector fields. We calculate 

\ the evolution of the level probability distribution as a function of time, and 

obtain rapid level diffusion over the energy levels. Our results imply a new 
mechanism of charmonium dissociation in QCD media. 
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I. INTRODUCTION. 



The propagation of the charged particles, both quantum and classical, through different 
types of media had long ago become an important branch of theoretical physics. (See e.g. 
an excellent review for the discussion on the present status of the subject). However, 
much less is known about the propagation of neutral bound systems and the wave packets 
with the zero total charge but nonzero higher multipole moments. This is especially true 
for quantum systems. Recently there was a revived interest in the subject. This is due 
to numerous possible applications in nuclear and particle physics, especially in the study 
of exclusive processes and the creation and the diagnostics of the quark-gluon plasma 0. 
Such systems can be also studied in electrodynamics, where it is possible to observe the 
propagation of e^e~ positronium through the random media. 

Although a significant advance in the study of the propagation of these objects was made 
recently, very little is known about the evolution of the fast wave the packets and the dipoles 
in the random media, and it's influence on the particle and nuclear cross-sections, diffractive 
processes and energy losses. 

The aim of this paper is to study the propagation of the highly energetic, but nonrela- 
tivistic (in it's cm. rest reference frame) electrically neutral dipoles and wave packets, with 
nonzero dipole moment, and with the internal interactions, through the random media. For 
simplicity, and in order to have an almost explicit solution we shall consider the simplest 
example of such kind — a harmonic oscillator moving through the media with the random 
vector fields. The latter can be any media, that where there are random vector fields, i.e. 
the fields transforming as vectors under the Lorentz group transformations. In nature these 
are electromagnetic and Yang-Mills fields. We shall call below any such media a vector 
media. The examples of the vector random media include amorphous solid state media [Q , 
usual electromagnetic plasma [Q, dense QCD media [^j , and quark-gluon plasma It 
is less clear what happens in the nuclear media or nuclei, where the fields are colorless and 
it is not obvious how they transform under Lorentz transformations. 

In our paper we consider the nonrelativistic oscillator, whose center of mass moves with 
the relativistic speed through the random media, which is characterized by the static random 
electric fields. In our simple model, the particles inside the dipole interact via harmonic 
oscillator potential (no Coulomb interactions between them), while the interaction with the 
external field is electromagnetic. However we shall argue below that at least qualitatively 
our results also hold for the interaction with the QCD colored fields, and for the general law 
of interaction between the constituents of the dipole. 

Our main result will be the development of the general formalism that will allow us to 
determine the density matrix, expansion rate and level diffusion for nonrelativistic harmonic 
oscillator moving with high velocity in the random media. We shall derive the effective action 
and the effective transition operator for the bound system propagation in the random media. 
This effective action will be used to study two closely connected problems: the propagation 
of the small initial size dipole/wave packet through the random media and of the oscillator in 
the ground state through the random media. We shall determine the energy level diffusion 
rate of the small wave packet for small and big times, as well as the probability that it 
occupies level n after time t. We shall also determine the expansion rate of the ground state 
and the probabilities that after time t the oscillator will stay in the ground state or will 
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move into the state with the energy level N. 

Taking into account that the electric fields transform as the components of the rank 2 
tensor during the Lorentz transformation, while the time as the component of the 4-vector 
we shall obtain for a certain interval of times the superdiffusion expansion law for the wave 
packet density radius: 

p' = {A,e + A,e)h' (1.1) 

Here ~ d — 1/d describes the quantum mechanical expansion rate of the wave packet 
andAi is the coefficient depending on the properties of the media, 7 is the Lorentz factor 
and,p is the radius of the dipole, and d is the inverse squared radius of the initial state 
(for the ground state d=l). Thus, random media may strongly infiuence the cross-sections 
and energy losses. The physical reason for "superdiffusion" is the extraction of the resonant 
(i.e. with oscillator proper frequency) mode from the spectrum of the random electric field 
fluctuations. 

This superdiffusion effect is very similar to the one discussed recently in refs. [^,^ in the 
context of the random walk of the classical particles in the time dependent random potential 
in statistical mechanics. 

Our results imply a new mechanism for charmonium dissociation in vector media:it is 
excited to higher energy levels due to the scattering on random electric fields, and then 



dissociates into DD pairs 10,11 



The calculation will be made in the eikonal approximation, i.e. we shall assume that the 
center of mass of the dipole moves along the straight line with the ultrarelativistic constant 
velocity Vq ~ c, where c is the speed of light. We shall neglect in this approximation both 
the deviations of the center of mass from the straight line, and the velocity changes due to 
the interaction with the random fields. The latter corrections will be taken into account 
elsewhere. We shall assume that the size of our dipole is small relative to the inhomogeneity 
scale of the media ^ 1/k. 

As it was mentioned above, all our calculations will be carried in the simplified model 
of the dipole with the oscillator interaction between the constituents and electromagnetic 
interaction with the external field. We neglect the Coulomb interaction between internal 
constituents of the dipole. We expect that our result will not change qualitatively if we 
generalize the external field to a color one, and take into account Coulomb interaction 
between the constituents. 

We shall see that the infiuence of the random media on the oscillator has very transparent 
mechanical analogue in the classical problem of the excitation of the oscillator by the external 



resonance force |T2|. The random media obviously has the role of this external force. On 
the other hand previous investigations, using this analogy, correspond to the study of the 
proper oscillations of the oscillator. 

Some of the aspects of the problem of the dipole propagation were considered before p , 
p!3| , p^ , but the latter authors did not take into account the infiuence of the media on the 



propagation of the dipole and it's expansion, and the change of it's quantum state during 
the propagation. They only considered the interaction free dipole, neglecting the infiuence 
of the media on the dipole structure and the internal interactions. Although as we shall 
see the internal interactions do not lead to new qualitative effects (as it was noted in ref. 
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p!3|), the interaction with the media does play a major role if time is not too small. The 



only previous attempt to take into account the influence of the media on the propagation of 
the oscillator in the media was made in ref. 0,0, but in the different context. The latter 
authors considered not vector, but scalar media, in particular nuclei, and parameterized it's 
influence by taking into account absorption. 

The problem of harmonic oscillator propagation naturally arises in the context of the 
propagation of the bound systems like charmonium through the quark-gluon plasma. 

Note that there are strictly speaking two different problems that can be discussed in the 
connection with the oscillator: one is the propagation of the two-particle bound state-dipole 
through the media, another is the propagation of the wave packet. In the approximation 
used in this paper we shall see that technically this is the same problem. The quantum 
dynamics of the bound state can be reduced to the quantum dynamics of the wave packet. 
So below, if not stated otherwise we shall not distinguish between them. However, beyond 
the approximation of the homogeneous electric field discussed below, these are two different 
problems, with different Lagrangians. 

The paper is organized as follows. In the second chapter we shall develop the general 
formalism for the description of the quantum dipole moving through the random media. 
First, we shall derive the effective Lagrangian for the propagating oscillator in this media, 
which we shall model by the usual plasma. We shall neglect the radiation. The discussion 
of this important phenomena will be considered in a separate paper. Next we shall consider 
the simplified model, that leads to the same qualitative results. We shall show that in 
the eikonal approximation it is sufficient to consider the behavior of the oscillator in the 
time-dependent homogeneous electric field. We shall calculate using the functional integral 
method the effective evolution operator and the probability to find an oscillator at particular 
levels as a function of time. In the third chapter we shall find the averaged density matrix 
and study the expansion of the small wave packets and see the infiuence of the media on the 
wave packet expansion. In the fourth chapter we shall consider the probability evolution for 
the ground state wave packets and apply our results to the simplest model of charmonium. 
We summarize our main results in conclusion. Some details of the calculations will be 
discussed in Appendices. 



II. THE MODEL OF THE OSCILLATOR EVOLUTION. 

A. The Formulation of the Problem. 

The most widespread model of the random media where the bound system can propagate 
is the conventional model of plasma 0] . Such model can be used qualitatively to describe 
the amorphous material. The media is represented by a random set of ions that generate 
electric potentials with the Debye screening. 

V{r) = g exp{—Kr)/r (2.1) 

Here g is the electric charge. If we neglect in this model three particle correlations, the 
potential— potential correlation function in such system is ( [Q) 
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L{r-7^) =< V{r)V{r) >= exp {- n\f - \) / \f - \ (2.2) 

A is the constant dependent on the properties of the media. We will assume that the 
correlations are gaussian, i.e. that all correlation functions can be expressed through the two- 
point correlation function ( p.2| ). This is the usual assumption in the statistical mechanics. 
More complicated systems with nongaussian correlations are usually highly nontractable. 

For the QCD media or the quark-gluon plasma the most wide-spread model is the model 
discussed in ref. |15|. This model is very similar to the conventional model of plasma, 



discussed above, except that instead of the conventional electric charges the particles have 
the color charges T"', that are the generators of the SU (3) color group, while the interactions 
contain factors g'^T°'T^ that are very similar to the charge products. If we take into account 
only the binary interactions, the plasma becomes color- diagonal, and we can neglect the 
color factors, that will not make qualitative influence on our results. The potentials and the 
correlation functions are the same as in eq. ( |2.2| ) (except some numerical factors whose 
influence will be discussed in a separate publication). 

Note that in terms of the language of the statistical mechanics for all cases (plasma, 
QCD media, quark-gluon plasma) our problem is the well known problem of the description 
of the propagation of the wave packet in the random gaussian media characterized by given 
correlation functions ( see Appendix A for a more detailed discussion). Consequently we 
consider in this paper all types of plasma as an electromagnetic plasma. 

Consider now the fast moving neutral system of two particles (dipole) with a harmonic 
interaction between it's components in the center of mass reference frame: 

W{\r-r\) = N\r-r'\^ (2.3) 

The system is moving in a random media discussed above with a velocity vq ~ c, where c 
is the speed of light. For this ultrarelativistic speed we can use the eikonal approximation. 
We assume that the center of mass moves along the straight line with a constant velocity. 
Radiation will be neglected. 

According to our assumption the oscillator is nonrelativistic, consequently we can sep- 
arate the motion of two constituents of the dipole into the motion of the center of mass 



and the relative motion of the two particles ( see ref. [0 and Appendix B for detailed 
discussion) : 

L = Lcm. + Ln.r. (2.4) 

In the eq. ( [2.4] ) the first term describes the free motion of the center of mass with the 
speed Vq, while the second is the Lagrangian of the relative nonrelativistic motion of the two 
components of the dipole in the cm. reference frame. The latter Lagrangian can be written 
for the case of the small dipole with the size <^ I/k, 

L„.,. = ^M{u^ - u^u') - g{V{n, t) - Vif2, t)) (2.5) 

Here V is the electric potential in the oscillator cm. frame, fi,f2 are the coordinates of the 
dipole components in this frame, and u = ri — r2 is the relative coordinate of two constituents 
of the dipole. The derivation of the Lagrangian (pT5|) is discussed in Appendix B. 
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For the small size dipole we may use a Taylor expansion: 

V{n,t) - V{f2,t) = (fi - rl)^^^^ (2.6) 

Here V(0,t) is the electric potential in the cm. reference frame at the moment t (t here 
is the laboratory time) that corresponds to the random potential in the point where the 
center of mass is located, it's derivative is obviously the electric field in this point. Thus for 
the small size dipole, with the size -C we can neglect inhomogeneity of the field and 
reduce the problem of the dipole motion to the problem of the oscillator in the homogeneous 
electric field 

(measured in the oscillator cm. frame in the point corresponding to the oscillator cm.). In 
this approximation, and this is an approximation that we will use in the present paper, the 
problem of the moving oscillator is equivalent to the problem of the propagation of the wave 
packet in the electric field E. For the larger dipole we must use the full second term in the 
Lagrangian ( |2.5D , taking into account the inhomogeneity of the field. 

The same separation can be made for the nonrelativistic neutral wave packet, whose 
center of mass moves with a relativistic speed. In this case the nonrelativistic Lagrangian 
that describes the space-time evolution of the dipole in it's rest frame can be written as 

Ln.r. = \M{r^ - uy) - gV{r, t) (2.8) 

Here r is the distance between two particles in the rest frame of the dipole or the radius 
vector in the cm. at rest reference frame for a wave packet. V is the potential in the wave 
packet cm. reference frame. The Lagrangian (p.5| ) in the approximation (|2.6| ) is evidently 
a particular case of eq. (|2.8|) . 

Let us go to the rest frame of the dipole. In the latter frame our problem is evidently a 
problem of the nonrelativistic oscillator with the internal interaction given by eq. (|2.3| and 
in the external electromagnetic field. 

Electromagnetic field acting on the particles that built up the oscillator can be easily 
found by making the Lorentz transformation |]T^. (Note that within the present approxi- 
mation, since the relative motion is nonrelativistic, the relative speeds of the two particles 
can be neglected to perform the Lorentz transformation according to cm. speed vq). 

After the Lorentz transformation (see appendix B for details), the laboratory frame 
4-potential (V, 0, 0, 0) of the electromagnetic fields becomes (K.m., 0, 0, A^): 

Km. = ^7, ^3 = -^7 

The electromagnetic tensor E,B [B = Q m the laboratory reference frame) becomes: 

-Ec.m.1,2 = -£'1,27 ) -^'c.m.S = -S'S 
^c.m.l = ^^0^27, ^c.m.2 = -^^0^l7, ^c.m.S = (2.9) 

Here 7 is the Lorentz factor: 
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7 = 1/v/l - ^11 (2.10) 

We see that the transverse electric field is greatly enhanced, while the longitudinal field 
remains the same. 

We also see that in addition to electric fields there arises a magnetic field equal to 

B = voxE (2.11) 

Fortunately, this magnetic field can be neglected. First, the magnetic field is of the 
same order as the electric one, but the Lorentz force acting on the dipole is suppressed as 
v/c « 1, where v is the relative speed of the dipole constituents, since we assume dipole to 
be a nonrelativistic system. Second, the magnetic field is the field with very short period and 
high frequency ~ K7. If we consider the classical motion in such periodic field, we shall come 
to the system of differential equations for the particle moving in combined oscillator-Lorentz 



force field, where we can use the Floque theorem. The Floque theorem |18| says that the 



system of the ordinary differential equations with periodic coefficients has the solution: 

x{t) = A{t)C{t) 

where B is the periodic matrix function and A is the diagonal matrix, with entries exp Xt , the 
coefficients A are called the eigenvalues of the problem. The matrixes are 6x6 and include 
both the 3 coordinates and its derivatives in time. Due to the energy conservation law and 
the fact that the classical energy of the oscillator is ~ {dx/dtY +u}qx'^, it is clear that all 
eigenvalues of the problem must be imaginary and appear in pairs (since the solution is real). 
Thus magnetic fields may contribute only to high frequency small fiuctuations along the 
classical trajectory corresponding to the motion without the magnetic field. Consequently, 
it's effect can be safely neglected. Clearly the same is true for the quantum mechanical case, 
since the corresponding path integral is saturated by the classical trajectories. 

We have seen that the fields that act on the oscillator in it's rest frame are time- 
dependent. In the eikonal approximation the coordinate of the bound state mass center 
is 

z = vot = VoTj, 

where T is the proper time in the oscillator frame, and 7 is a Lorentz factor: 

T = t/7 (2.12) 

Using the time and frequency transformation and field transformation laws we can write the 
correlation function of the potentials in the oscillator rest frame. The correlation function 
will become time-dependent: 

< V{f, T)V{i^, T') >= A^-f^ exp{-Kq)/q (2.13) 



q = ^{x - x'Y + {y- y'Y + 72(2 -z' + vo{T - T'))^ (2.14) 

In the latter equations x and y are the transverse coordinates that are the same in both 
reference frames, while z-s are the z-coordinates of the oscillator in the cm. frame. Due to 
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the Lorentz contraction, we can neglect z dependence and therefore assume z = z' (this can 
be done in Quantum Mechanics, where we neglect radiation). Recall that T and T' are the 
proper time of the oscillator. Thus we see that our problem of the moving oscillator wave 
packet has been reduced to the problem of the nonrelativistic oscillator at rest in the high 
frequency external electric field that is strong in the transverse plane, and it's longitudinal 
component has the same strength as in the media at rest, with the random field correlator 

In the small oscillator limit discussed above, the electric field can be considered as ho- 
mogeneous, and we get the correlation function 

< V{J)V{J') >= AVexp(-Kg)/g (2.15) 



q = Vo-f\T -T'\^-f\T -T'\ (2.16) 

Here and below except in the 7 factors we assume Vq = c, and use the units where c = 1. 

We shall argue below that the results obtained in the homogeneous field approximation 
hold qualitatively even for the large scale oscillator. 

The eq. ( p.l4|) holds also for the dipole, since in this case we work from the beginning 
in the homogeneous field approximation (|2.7| ). 



B. The General Formalism. 

Now we have to consider the problem of the nonrelativistic oscillator in the random 
media, with the correlation function of the scalar potentials defined by eq. ( |2.15| ). This 
problem is very similar to the problem of the propagation of the electromagnetic wave in 
the random media characterized by the gaussian distribution function (see e.g. ref. ||19||). 

Our main goal is to find the level distribution of the oscillator at the moment of time 
t if at the moment he was in the state characterized by the wave function \Ef(r, 0), or 
equivalently by a set of coefficients a„(0): 

^(f,0) = Sa„^„(r) (2.17) 

Here are the wave functions of the stationary states of the pure oscillator without the 
random fields: 



^„(x) = (Ma;o/vr)^/'y^exp(-Mu;oxV2) (2.18) 

(we shall choose the units where the Planck constant h = 1). The probability that at the 
time t the oscillator will be at the level n is given by 

P(n,t) =< |a„(t)|2> (2.19) 

The average is over all field configurations. The simplest way to calculate the latter average 
for the oscillator is to calculate the average product of the Green functions: 



7 



F{xo,yo,xi,yi,t) =< G{xo,xi,t)G*{yo,yi,t) > 



(2.20) 



Here G is the oscillator Green function in the external field V, that describes time evolution 
of the oscillator in the given realization of the field V. The average is taken over all external 
fields V with a Gaussian weight. Once we know the function G{xo,yo] Xi,yi,t), the time 
evolution of the a„ coefficients can be easily calculated: 

\an{t)\^ = J d^xod^yod^x,d^yiF{xo,yo,xuyi,t)^{x)^{y)^n{x')^n{y') (2.21) 

Here we assumed that at t=o the oscillator was in the state described by the wave function 
^(x,0). 

One possible way to carry the calculations is to use the functional integrals method. The 
Green function of the 3D quantum oscillator is well known and can be represented as (see 
e.g. ref. EO], PTI, p2| and references therein): 



G{x,y,t) = f dx{t)exp{i f ]-M{x^ - ujlx^))ds (2.22) 
J Jo 2 

Here the integration is over all paths with 

x(t) = xi,x(0) = xo (2.23) 
Consider now the arbitrary wave packet. The full functional integral we are interested in is 

F{xQ,yo,xi,yi,t) = j dx{t) j dy{t) j dV{x,t)e-x^{i{^{{xf - ujIx"^) -V{x) 
-{^^{{yf-u;y)-V{y))) 

X exp(— / / dsds'd^xd^yV(x, s)K(x, s;y, s')V(y, s')) 
Jo Jo 

(2.24) 

The operator K is just the inverse of the correlation function of the potentials V, and is 
determined by the properties of the media. The integration is over all paths with 

x{t) = xuy{t) = yi; x{0) = xo, y{0) = yi (2.25) 



Note that in the eq. ( |2.24| ) the integration over the oscillator paths is in Minkowski space, 
with the usual Minkowski time t, and i before the action. However, in the part of the 
functional integral connected with the integration over V there is no i before the action, and 
the integration is purely statistical one. This integration is analogous to Euclidean space 
functional integral in the Quantum Field Theory, and it is evident that it has the same 
nature as the coupling of oscillator to the finite temperature thermal bath. 

It is very hard to carry out the exact functional integration in eq. ( |2.24| ) even though 
the corresponding functional integrals over V and x,y are Gaussian. Indeed, integration 
over V can be carried out. We can proceed in the standard way, as in the theory of the 
random fields, by looking for extremum V trajectory, and then substituting the result into 
the functional integral ( p^.24D that will be now the integral only over x,y. In order to carry 
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the integration, we write the term V(x(t),t)) in the action as / xS{x — x{t))V {x , t). Varying 
the action over V(x,t), we obtain: 

i6{x - x{t) - i6{y - y{t)) = - /* K{x, t; x't')V{x', t') (2.26) 

Jo 

This is the equation for the extremal trajectory. The latter equation (except the factor i 
due to the fact that we carry the average in the Minkowski space) looks like the potential 
created by two opposite charges moving over the trajectories x(t) and y(t). The equation 
can be inverted giving 

V{x,t) = -i ds{R{x,t;x{s),s) - R{x,t,y{s),s)/2 (2.27) 



Here R is the correlation function (|2.14|) , that is inverse to the operator K. Substituting this 



expression for V in the functional integral ( p^.24D we get 



F{x,y,x',y',t)= dx{t) dy{t) expt{M/2){x' - y' - uox' + Uoy') - SA{x{t),y{t)) (2.28) 



Here the effective action A is 

Sa= f dt' f ds{R{x{t'),t';x{s),s) + R{y{t'),t';y{s),s) 
Jo Jo 

- R{x{t'),t'- y{s), s) - R{y{t')- 1'; xi-s), s)) (2.29) 

Note that the above action looks quite similar to the famous Feynman action for the polaron 
p3| . It describes two particles with the singular interaction potential, that includes both 
interaction with it's own trajectory and the interaction with the trajectory of the second 
particle-the situation quite similar to the problem of electrostatic dipole. 

It is possible to deal with the latter action using mean field approximation or perturbation 
theory, however for fast particle the problem can be significantly simplified. 

In order to simplify the problem further, first note that the main contribution into the 
action comes from the region t' ^ s, due to the exponential Yukawa cutoff in G. In the area 
t' ^ s, we can expand 

f(t') - x{s) ^ -{dx/dt'){t' - s). 
Then the radius in the Yukawa correlation function |2.14| will have the form 



^{f - sY-i^{vl + {dx/dt'Y) (2.30) 



Since the system is nonrelativistic, the second term in the brackets in eq. (|2.3CI| ) can be 



neglected, and the system reduces to the usual quantum mechanics, i.e. we can neglect 
both the longitudinal and the transverse self-motion. Concerning the interaction with the 
trajectory of the second particle, the contribution is clearly maximal when x(t) is close to 
y(t), and in this case also we can neglect the transverse motion. This means that for a 
relativistic particle in the first approximation we can neglect the space dependence of the 
electric field and assume that the system is in the homogeneous time dependent random 
electric field. In fact it was clear from the very beginning that we can take our electric 
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field as homogeneous. Indeed, if we look at the action of the electric field, we see that 
the term due to the space inhomogeneity contributes into the action as k^E"^, relative to 
/£2^2^2 time evolution, and thus can be neglected, even if the initial wave packet 

is not of the small radius. We shall see below, that the inhomogeneity of the electric 
field must be taken into account only at the times of the order of the minimum of two 
scales: 'j^^^/ {uquY^^, 7^/ \J{d — l/dhzuo, when the longitudinal expansion of the wave packet 
becomes significant— of the order of I/k. Thus our results can be applied even for the more 
general case when 1/k < I/uq. 

Concerning the small dipole, we deal with it in the approximation of the homogeneous 
electric field from the beginning, the field being equal to the field in the cm. of the dipole. 

Next, we need to deal with the Coulomb singularity of the Yukawa potential. We shall 
follow the line of ref. The contribution of the singular part of the interaction to the 

functional integral is zero: the probability of such configuration is zero. Moreover, it is 
known that such singularity really does not appear due to quantum fluctuations, and it's 
appearance is due to the quasiclassical method of the path integral calculation, since in 
the classical case there can be fall on the center. Consequently, if we are interested in the 
qualitative results, we can substitute the singular correlation function with a nonsingular one, 
omitting the coulomb part and dealing directly with the electric fields. As for the strength 
of the electric field we can take the strength of the field it is clear from eqs. ( |2.5| ),( p77D that 



we have to take the field in the reference frame where the center of mass of the system is at 
rest. 

Thus we come to the conclusion that both qualitatively in the leading eikonal approxi- 
mation and quantatively we can describe our model by an oscillator coupled to the homoge- 
neous time dependent random electric field with the correlation function in the longitudinal 
direction given by 

< E,{t)E,{s) >= exp{-Kjvo\T - s\) (2.31) 

and in the transverse plane with the correlation function 

< Eiit)Ejis) >= (527^) exp(-fi:7i;o|T - s\)5ij (2.32) 

Here T and s denote the proper time (in the moving oscillator frame of reference), and 
the frequency is increased greatly by 7-a Lorentz factor. The corresponding potential is 
V{x) = Ex, and the corresponding Lagrange density in the action can be rewritten as 

SSa = A^j dtE(t) {d^ + K^-i^)E{t)dt (2.33) 

Note that the simplifications we did are very similar to the Feynman trick for polaron, where 
he substituted the complicated nonlinear action of a polaron by an action of the harmonic 
oscillator. 



C. The Calculation of the Effective Action. 

Now we shall calculate the 4-point function F in the above model. In order to average 
over the field configurations we shall use the action 
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S = A^j {{dEf + K^^^E^)dt (2.34) 

The constant A is determined using the relation with the correlation function of the oscilla- 
tor: 

< E{t)E{t') >= B\l, 7^) exp -(|t - t'\K-f) (2.35) 

(the brackets depend on whether the field is in longitudinal or transverse direction). Then 
for the longitudinal field we have : 

EliO) =Bl = l/iAli2K^)) (2.36) 

For the transverse field we have ||2l| : 



E^iO)=B^Y = l/iAii2K^)) (2.37) 

Here both electric fields El and E^ are the fields in the frame moving with the velocity 
f ~ c where the harmonic oscillator is at rest , so that B is the \/El{0), where index L 
means the laboratory system, and the time t here and below refers to the frame moving 
with the oscillator. Hence for the longitudinal field 



For the transverse field: 



AL = l/iBLi^'/^^V2) (2.38) 



At = 1/{Bt1^^'^^V2) (2.39) 



Let us now repeat the derivation of the 4-point function G for the concrete example of 
the averaging weight E'^. The derivation is slightly different from the previous chapter, since 
there we dealt with potential V and varied over it, while for the final simplified model we 
can deal directly with electric field and vary over it that simplifies the problem. 

Due to the obvious factorization we can start from the one- dimensional problem. The 
whole calculation is carried in the rest frame of the moving oscillator. 

More explicitly, the functional integral we need to calculate is 

F{x, y, x', y', T) = J dx{t)dy{t)dE{t) exp iS{x, ,y,E)- S{E) (2.40) 
Here S(x,y,E) is the action 

S = / dt— (i^ - uj^x^ -y^ + uW) -Ex + Ey (2.41) 
Jo 2 



The statistical weight S(E) is given by eq. (|2.34|) . The boundary conditions for x,y are 



evident. Note that the first term in the functional integral ( |3.9| ) is multiplied by i while the 
second is not. We first will carry the integration over the random fields, and then consider 
the effective action for x and y. The integral over the electric field E is Gaussian and can 
be taken. We first need to vary the action S over the fields E and then substitute the result 
into the action. Varying over E we get: 
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i{-Xi{t) + y,{t)) = A'id" - K'Y)Ei{t) (2.42) 
The latter equation can be immediately solved using the Green function of the operator 

(a^ - K^) exp{-\Kt - t'\)/{2K) = 6{t - t') (2.43) 
Consequently, we write 

Ei{t) = -(i/(2A2«7)) rexp)(-|i - s\){xi{s) - y,{s))ds (2.44) 

^0 



We now substitute the above expression for E into the action in order to get effective 
action for x and y. We have for the terms that previously contained the electric field 

82^- r r dtdsex.p(-K,-f\t - s\)(x(t) - x(s)(y(t) - y(s))/(2A'^K-f) (2.45) 
Jo Jo 

Note that the internal integration over s in the action goes up to T and not t, like in the 
Minkowski field case. The reason is that the averaging over field configurations takes part 
only after the concrete trajectory in x(t) is realized. Physically, we first calculate functional 
integral for given E reahzation, and then average over all E. 
We now have the effective action -^i + 5'2, where 

Si = iM{x^ -f- ujqx^ + ujoy'^) /2 (2.46) 

The calculation of the functional integral proceeds in the standard way. We seek the solution 
of the classical equation of motion, giving the exponent, and then calculate the preexponen- 
tial factor. 

Let us start from the equations of motion. We get 

iM(-^ + uj^x) = - ^ exp(-K7|t - s\ {x{s) - y(s))/(A'«7) (2.47) 



iM{-^+uoy) = - ^ exp(-/^7|i - s\){x{s) - y{s)/{A''K^) (2.48) 
We see that the equations for x and y are identical, thus 

^-^^ + u;lix-y) = (2.49) 

The boundary conditions are: 

Si = x{T) y{T) = xi- yi- 60 = x{0) - y{0) = xo - yo (2.50) 
The corresponding solution is evidently 

x{t) - y{t) = {Sosm{uJo{T - t)) + 5i sin uot)/ {sin uT) (2.51) 
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Once we know x— y, we can substitute it to eq. ( |2.48|) and get the equation for y, which is 
just the equation for the oscillator under the action of the external imaginary force. 

^ + ujIx = i/{{MA\-f) sin(cJoT)(fi:V + ^o)) 

X (exp(— K7t)(co'o5i — So{uJo cos{ljqT) + K7sin(co'oT)) 

+ exp(-K7(T - t)){5oUJo - 6i{uJo cos{uJoT) + ^7 sin(u;oT)) + 

+ 2k7(5i sin(cjot) - 5o sin(cJo(t - T))) (2.52) 

Note that the external force that appears in the r.h.s. of the eq. ( |2.52| ) has three parts. 
Two of them are exponentially suppressed and are different from zero only close to t=0 or 
t=T. However, the third term is not suppressed, and, moreover, is a resonance force applied 
to the harmonic oscillator with a frequency ujq. This is of course what was to be expected, 
since the random field distribution contains all possible frequencies, including the resonant 
one. The solution for x for given boundary conditions will be a sum of the solution of the 
homogeneous equation plus the solution of the nonhomogeneous one: 



1 /•* 

x(t) = — / sin(ci;o(s — t))u{s)ds + Xocos{ujot) + Bsin^uot) (2.53) 
ujq Jo 



Here u(s) is the r.h.s. of eq. (|2.52| ). The coefficient B(T) is determined from the boundary 
conditions. 

The general solution of the nonhomogeneous equation is a sum of five terms: 

X{t) = Xa{t) + Xb{t) + X,{t) + Xd{t) + Xe{t) (2.54) 

The first term is due to the resonant part of the external field and is given by 

Xa(t) = iL{6iK'yt cos{u!ot) — 6oK'yt cos{ujq(T — t)))/u!o (2.55) 

The second part includes the terms that come from the integration of the resonant part that 
are not linearly enhanced as the term above: 

Xb{t) = 'iL—{5Q sin(c<Jot) cos(co'oT) — 5i sin(u;o'^)) (2.56) 

Here 

L = l/{A^K-i)Msin{tooT){K^-i'^ + uD) 

The third part comes from the integration of the nonresonant force, but is not exponentially 
suppressed, however it is suppressed by additional power of K7 relative to the previous term: 

Xc{t) = (^-^/(^7^ + ^o)){^o^i — 5o{ujq cos{uJoT) + ^7 sin(u;or))(u; cos(u;o^) — K'-^ sm.{ujQt)) 

(2.57) 

Finally the 4th term is the exponentially suppressed part coming from the nonresonant part 
of the external force: 
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Xd{t) = (iL/(K7^ + ujl)){-iUo{iUo5i - 5o{uJo cos{uJoT) + ^7 sin(u;oT)) exp(-K7t) 

( exp(— k7T)(k7 sin(c<Jot) + cuo cos(c<jt)) — c<j exp(— ((T — t)))) (2.58) 
Finally there is a solution of the homogeneous equation that can be rewritten as: 

Xe{t) = Xo COs{uJot) + (Xi - Xa{T) - Xb{T) - Xc{T) - Xd{T) 

— Xe(T) — XoCOs{ljqT)) sin(a;o^)/ sin(a;o^) (2.59) 
Note that if we are interested in time scales 

tL > 1/k 

we can neglect the terms that are exponentially suppressed. The terms in x^ have the same 
structure essentially as in Xb{t) but are enhanced by at least K'-f/ujQ. Hence we can write for 
tL > 1/k 

X{t) = Xa{t) + Xbit) + Xe{t) (2.60) 

Here in Xe{t) we can in the expression for B(T) neglect all the terms except coming from 
Xa(T) + XbiT) + xq cos^ujqT). Then for B(T) we have 

B(T) = {xi — a;o cos(co'oT) — iLK'y(T{6i cos(co'oT) — 60) + sm{LjT){6oCOs{LjQT) — 5i)/ujq)/ sin(a;o^) 

(2.61) 

We see that the classical trajectory of the oscillator in random field corresponds to the sum of 
two terms: first, the usual harmonic oscillations (due to the homogenic part of the solution), 
second the linear expansion (although suppressed as I/7) in the imaginary direction . 

We now calculate the action as the function of the boundary conditions using the solution 
( p.60| ). Due to the equations of motion the action is just 

i{x{T)x^ - y{T)y^ - Xox(O) + mVo), (2.62) 

The other terms in the action disappear due to the equations of motion. It is easy to obtain: 

S = ^^{{xl -yl-xl- yl) cot{ujQT)) - 2{xqXi - y^yi)/ sin(u;oT))) 
+ i{LK-i / sm{ujQT)){{5l + 5l){ujT - sin(2tuT)/2) + 

+ 25o5i(sin(co'o7') — ujqT cos{ujqT)) (2.63) 

This action gives us the transition operator as a function of boundary conditions. Note 
that it contains three types of terms: 1) real part; 2) Imaginary part with linearly enhanced 
terms due to quantum fluctuations 3) Imaginary part with only trigonometric terms. Note 
that we neglected for all the terms that are suppressed by exp(— T)) and by the powers of 
7. The transition operator in the quasiclassical approximation will be 

F = exp{iS) (2.64) 
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We need now to consider the preexponential factor, and normalize to the evolution oper- 
ator of the harmonic oscillator without electric fields. In order to calculate the preexponent 
we need to expand: 



x{t) = Xci{t) + ^ c„ sin(27rnt/T) (2.65) 





y{t) = Vciit) +J2bn sin(27rnt/T) (2.66) 



Since the cross terms with Xd evidently cancel we need to consider the integral in the action 
due to the sum over n. It is equal to: 

P = ^exp(^M(6^ - 4)((27rn/T)2 + c.^) 

I rT rT 

- ^JI^ E(cn - bn){cm - sm{27rns/T) 

X exp(— K7|t — s\ sm{2nmt/T)dtds (2.67) 

Wc can now move to the integration over the variables Un = bn — Cn and Vn = bn + Cn- Then 
we evidently get after the integration over Vn the product of delta-functions 

l[S{ur,{{2nn/Tf + u;l)) 

Then we can evidently let Un = 0- The preexponential factor P will be the same as the 
corresponding factor for the oscillator without the electric field which is evidently: 

n=oo 

Now we can calculate the transition operator G, whose matrix elements give the transition 
probabilities we are looking for. We obtain 

F{xo,yo;xr,yi;T) = . / cxp((iMa;/2)( . ^ ((xg + xj - yl - yl) 
cos {ujoT) - 2{xoXi - yoyi))) 

X exp -{MLK-i . } A b\ + bl){u)^T - sin(2a;oT)/2) 
sm(^c<^o-' ) 

+ 25o(5i(sin(a;or) -u;o7'cos(u;or))) (2.69) 

Note that the first term is just the usual evolution operator for the harmonic oscillator, while 
the second term contains the action of the random field. The operator is written in terms of 
the proper time, in the oscillator reference frame. In order to work in the laboratory frame 
we need to substitute T = tL/7, where ti^ is the laboratory time. We can use the latter 
transition operator to calculate transition amplitudes as the functions of time. 
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D. Transition Amplitudes. 

Once we know the transition operator ( |2.69| ) we can calculate the evolution of the initial 
distribution. If in the initial state t=0 the oscillator is in the state k, the probability that 
after time t the oscillator will be in the state k is 

|a^(T)p = J dxidxodyidyoHnixi)Hn{yi)Hkixo)Hk{yo) 

exp -{Mujo/2){xl + yl + xi + yl))F{xo, yo, x,, yi, T) (2.70) 

Here iV„ is the normalization factor: 

iV„ = (Ma;o/vr)'/VV(2"^!) (2-71) 

In order to carry the actual calculation we need to go to the variables Sq = Xo + yo, Si = 
^1 + yi] i^O) ^1 and carry out the corresponding gaussian integral over the four variables. Note 
that due to the cross-term between 6o6i in the transition operator F we have 

|<a>|V<l«P> (2.72) 

The easiest example to calculate is the transition probability |aQp for the oscillator en- 
tering the media in the initial state, and then staying in the initial state. The corresponding 
gaussian integral is 

\aUT)\' = l/VH (2.73) 

where H is given by 

H = l + 2uoDtL + AD\ultl - 7^ sin\uoQUh)) (2.74) 

Here 

Dl = 2BlK/{Muol{K''-i^ + uD) (2.75) 
for the longitudinal motion and 

Dt = 2B^K-fy{Mujl){K^j^ + cuD) (2.76) 
for the transverse motion. We see that for >> •y/u the dependence is 

^ 1/iDuJotL) 

The result for the arbitrary n,k is given by the derivative: 

< |a^(T)|2 > = H,,((dx, + dx,)/2) X R^Hdx, - 9aJ/2) x HkUd,, - dfi,)/2) 

Hk{d^, + d^,)/2)I{T,fi, = 0,/i2 = 0, Ai = OA2 = 0) (2.77) 

Here the integral I is 
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J dsodsid5od5i exp{iS + hiSq + H25o + Xi5o + \25i - MuJo{5f + 5^ + sl + si) /A)/ sin(co'oT) 

(2.78) 



The integral I is a gaussian integral that can be easily taken: 

1 



Here 



exp{-Z/H - {Xi + fii)/MuJo))) (2.79) 



Z = {Muj{l + D{ujT - sin(2cjT)/2))(/i^ + A^) 

- (Mcu(l + D{ujT + sin(2cuT)/2))(/i^ + A^) 
+ 2iL)Mu;sin2(cuT)(/ii/i2 + A1A2) 

- 2i(A2/ii + Ai/i2) sin(cuT)Mcj(l + DujT) + 2/iiAiMcu(cos(wT) + DujT + D sin(cjT)) 
+ /i2A2Mtu(cos(wT) + DujT - D sm{ujT)))/2H (2.80) 

The problem is of course SDimensional, and the explicit expression for 3D coefficients 
are just a direct product of one-dimensional coefficients in x — y plane and in z direction 

= n <c (2-81) 

1=1,2 

Here we took into account that for sufficiently big 7 the transverse fields are strongly 
enhanced and we need only to take into account the transverse transitions. 

The evolution of the arbitrary initial distribution ( p.l7|) can be calculated in the same 
way. 



III. EVOLUTION OF THE ARBITRARY WAVE PACKETS. 

In the previous section we considered the time evolution of the arbitrary initial states. 
In this section we shall consider the evolution of the small wave packet that has at t=0 the 
form in the transverse plane: 

^{y) = {2dMuJo/7ry/^ exp{-dMuJoy^) (3.1) 

We first calculate explicitly the density matrix 

p{x,y) =< ^{x,ty^{y,t) > 

where the average means the averaging over all random field configurations, and then use 
this expression to calculate the probabilities. Substituting eq. (|0| ) into the equations of 
the previous section, we get for the density matrix 

p{x,y,t) = exp{-6^{d/lQ + D^d{iult^ - -f^ sin^iuotf-f) + {ut + sm{2LUot/-f) /2.)D /A 
+ d'^D{uJot - sin(2tUot/7)/2) - ds^lG + 2iSs{{d'^ - 1) sin(2cJot/7)/32 

+ Ndj sin\uj,th)lA))/D{t))/^D{t) (3.2) 
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Here D is given by eq. ( p.76| ) and the function D{t) is given by 



D{t) = {d^ sin(cJot/7)^ + cos(cuot/7)) 74. + Dd{uJot - -fsin{2uJot/-f) /2) (3.3) 

We can determine density radius of the dipole/wave packet. The density radius char- 
acterizes the decrease of < > as a function of x^. Note, that this is completely 
different physical quantity than the effective radius of the dipole/wave packet, that deter- 
mines the effective cross- sect ions. 

The expansion of the wave packet density is determined by p{x,x) =< ^(s, t)*\E'(a;, t) > 

^ D{t)/d = (rfsin(cJot/7)^ + {l/d) cos(tuot/7))/4. + D{uJot - 7sm(tuot2/7)/2) (3.4) 
Note that for the small t<^ 'j/uj we have 

x^t) ^ {l/d +lid- l/d)uJoty^^ + 16Du;lty^^)/4: (3.5) 

We see that there are three pieces in the latter equation: first, the initial radius of the 
wave packet, second-the geometrical factor, that is equal to zero when the particle is in the 
ground state, and is bigger, the smaller is the wave packet radius (in correspondence with 
the uncertainty principle). This factor describes the wave-packet diffusion, that is present 
when there are no interactions. Finally, the third term is wave-packet radius independent 
and describes the density expansion due to a strong field. We see that the latter for small t 
goes as t^, and for the case of the ground state d = 1 this is the only term that contributes 
to the expansion. 

Before the time reaches tk ^ {d — l/d)/{uJoD) the density radius squared increases 
according to a square law, while for times > 7/u;o it increases according to the linear law. 
We can determine, till what time scale we can consider the field E as homogeneous. 



It is well known in the theory of random systems that the field can be considered 
homogeneous if the inhomogeneity scale in transverse direction is much bigger than time 
( and z) correlation scale. In our case the latter scales are 1/k and 1/(k7) correspondingly. 
However, due to the longitudinal diffusion the longitudinal-time correlation scale increases. 
Due to the fact that the longitudinal electric field does not change strength due to Lorentz 
boost, the expansion for the small packet in the longitudinal direction is determined by the 
geometrical factor. For the big packet with the radius of the ground state or bigger the 
expansion rate will still be determined by the diffusion. The longitudinal scale will reach 
the radius 1/k at the minimum of the times 7^/ \J{.d — 1/ d)ujQK, •y'^^^ / {u^kY^^ . This means 
that for all realistic values of d and for sufficiently large 7, our results, derived for the 
homogeneous fields, certainly hold up to the times of ~ 7/u;, that is the characteristic time 
in the interaction of the dipole and the target. 

For the case of the oscillator the same result can be obtained in the different way, using 
the Erenfest theorem. Indeed 

^ + u;^x = E{t)/M (3.6) 
This operator equation is easily solved: 
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x{t) = Acosiuti/^) + f dssiniuoit- s))E{s)'/M (3.7) 

Jo 

Then 

< p^{t) >= {B^K-f)) r dsds'sm{uj{t-s))smuj{t-s') < E,{s)E,{s') > (3.8) 
JO Jo 

for longitudinal diffusion, and 

< x^{t) >= (B^-f^K-f) / / dsds' sin(u;(t — s)) smuj(t — s') < Et{s)Et{s) > (3.9) 

Jo Jo 

for transverse diffusion. In this way we obtain the law of expansion: 

< p^(t) >^ B^{{2tL - 7 sin(2cut/7)/(2cj))) x /tV/(«;V + ^^) (3.10) 

We depicted the characteristic dependence of the square radius of the oscillator in the 
initial ground state and in the initial small wave packet state in figures ||,PJl^ and p!l| , |T^ ,|T3 
respectively. Note that in this derivation it is irrelevant if we consider classical or quantum 
problem. The coincidence of the results of two approaches is of course due to the Erenfest 
theorem. 

Once we know the laws of quantum expansion we can calculate the probability distribu- 
tion for the system to be in the n-th state. The probability is given by 

< |a„(t)p >= S{tr/yR{tY'^+'^/'P„Mt)/i2Dit)^S{t)Rit)) (3.11) 

Here 

Q{t) = D^d{ujH^ - 7^ sm{ujt/-ff) + D{ujt + 7 sin(2u;t/7)/2)/4 + 

N d'^{{ujt-sm{2ujt/-f)/2)/A. (3.12) 



Pn are the Legendre Polynomials, and D(t) is given by eq. (|3.3|) ; 

S{t) = -1/4 - {D^d'^iooH^ - sin2(cjt/7) cos\ujt/-f) + d^lQ 

+ Dd{ujt + 7 sin(4 * cjt/7) /4) /4 + Dd^{ujt - sin(4cjt/7) /4.)/4 
+ {d^ - l)2sin(2a;t/7)V64)/(4D(t)2)) 

( d/8 + D'^d{ujlt^ - 7^ sin2(cut/7)) + D(cjt + 7 sin(2c<;t/7)/2)/4 
D d^{ujt - sin(2cjt/7) /2) /4.)/D(t) (3.13) 



R[t) = +1/4 + {D'^d^iuoh'^ - sin2(cjt/7) cos^{uot/-i) + d^/lQ 

+ Dd{ujt + 7 sin(4 * ojth) /4) /4 + Dd^{ujt - sin(4cjt/7) /4.) /4 
+ {d^ - lfsm{2uthf/U)/{W{tf)) 

( rf/8 + D'^d{ujlt^ - 7^ sin2(cut/7)) + L)(cjt + 7 sin(2cut/7)/2)/4 
D d^{ujt - sin(2cjt/7)/2)/4.)/L)(t) (3.14) 

The sample graphs are given below ( see figures 00,0). 
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IV. THE GROUND STATE EVOLUTION. 



Consider now the oscillator entering the random media in the ground state. Then we 
can put d=l in the formulae of the preceding chapter and all results will be significantly 
simplified: 

R{t) = H{t)/{AD{t)), S{t) = D\ujie - f sm{uothf/D{t) (4.1) 

The probability is given by 

< M)? > = ((2Z^)^ K¥ - f sin(^ot/7)^)"/V(/ /(t))("+^)/^ 

X PniD^iooit^ - 7^ sinKt/7) V^(i)) (4-2) 

The sample graphs are given below (see figs. 

These results can be applied to the simplest model of charmonium moving through the 
Quark-Gluon plasma. This picture of charmonium is of course not realistic , since it does 
not contain absorption P, p!4| , p5| , coulomb interactions |TO| , p!T| , color and radiation. It is 



given here for illustrative purposes only. The more realistic picture of charmonium will 
be discussed elsewhere. In our model the charmonium enters the media (say quark-gluon 
plasma or nuclear matter) as the dipole in the ground state. The charmonium is scattered 
by the random electric fields in the media and is excited to the states with = 2 or higher 
states, where it decays via DD pair decay mechanism. We assume for simplicity that the 
decay probability is one, i.e. once the charmonium crosses the DD threshold [|TU|,|TT|, the 



total probability that charmonium will still be charmonium and not decay will be given by 
the sum 

i+j+k<No 

p{T)= y: i«(o-^jfc)r (4.3) 

i=0,j=0,fe=0 

Here A'^o is the energy level corresponding to the DD threshold. It is enough to take into 
account only transitions in transverse directions, so the sum is given by 

P{t)= E \«\"\ (4-4) 

p+s<No 

The result is depicted in figure |^. In the calculations we have taken the realistic ujq 
= 1 GeV. For the realistic k we take 0.12 GeV, for 7 ^ 30 (that corresponds to quick 
charmonium produced in the Fermilab experiment), and for the field B we take the field in 
the QCD vacuum |^ , that is say 0.11 GeV^, and for illustrative purposes we took Nq = 2. 

Our results imply that this mechanism can make a significant contribution in the disso- 
ciation of the rapid charmonium in the QCD media. 

We also did not take into account that there exist, for the small wave packet, backward 
transitions from higher levels to low ones, that must be discarded in the realistic model 
(since charmonium will dissociate). This question will be dealt separately. 
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V. CONCLUSION. 



We have studied the quahtative influence of the random media on the relativistic prop- 
agation of the nonrelativistic dipoles and wave packets with internal interaction ( see eqs. 



?2D, (|3.5|) ). We showed that the media: 1) enhances the expansion rate and leads to the 
new regime of the superdiffusion expansion at the rate ~ for a ground state and for 
the small wave packet for a sufficiently large time intervals if a field is sufficiently large; 2) 
leads to the diffusion over levels, leading to new types of diffractive events. 

We have seen that the wave packet expansion rate depends on the field strength and can 
be used to study the properties of the media. Moreover, our results strongly suggest that the 
energy losses and radiation of the dipole or wave packet propagating through the random 
media will strongly depend on the media (compare with the corresponding calculation for 
a color parton |2^), and can be used for the diagnostics of the quark-gluon plasma. Our 
results also imply an additional mechanism of the charmonium output suppression relative 
to the known ones [p5| , p7| . 

Our results imply the possibility of an additional mechanism of dissociation of rapid 
charmonium by excitation to the higher levels due to scattering on the random fields, and 
consequent dissociation into DD pairs. 

We have calculated the transition operator ( p^.69| ),( prUBD ) for the oscillator, and cal- 
culated the time dependence of the probabilities of the oscillator occupying the 
n-th level if it was in some initial state. The corresponding formula is given by eqs. 
(CT),(F^,(CT ) ,(|33^ ) ,(ro)- We used these results to study the simple illustrative model 



of charmonium, and were able to determine how many charmoniums will stay after time t 
(or, equivalently, if a charmonium beam goes through the shell of the width L = Vot) (see 
eqs. (0), (H). 

The results obtained in this paper were derived for a particular case of the moving 
oscillator. We expect that the same results still hold, at least qualitatively, for more general 
system, such as a Coulomb one. and. Indeed, the mechanism of superdiffusion expansion 
is the pick-up by the system of the resonance frequencies from the media field frequency 
spectrum. In the system other than oscillator, say in the Coulomb system, it will pick 
the frequencies corresponding to the distances between nearby levels, and consequently will 
quickly expand, as in the simple oscillator case. It will be very interesting and necessary for 
the study of realistic systems to take into account absorption P,p|. 



We have calculated the effective action and transition operator (|2.45| ) for the propagation 
of the general relativistic dipole in the media. This action, obtained after integration out of 
the electric field is quite general, and can be used for the arbitrary bound systems. However, 
for the cases other than harmonic oscillator it is hard to carry the exact calculations and 
one may resort to numerical methods. 

It will be very interesting to use our method to study a number of problems not considered 
here: the propagation of positronium through amorphous media and it's ionization rate, and 
the more realistic charmonium model that includes color and absorption, as well as to the 
detailed analysis of the infiuence of the nonhomogeneity. It will be also very interesting to 
carry the calculation of the radiation loss and it's infiuence on the dynamics of the dipole. 
At this point, due to the appearance of the traces of the several color matrices, we shall 
expect real technical difference between calculations in the usual plasma and quark gluon 
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plasma. 
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APPENDIX A: THE PLASMA MODEL AND IT'S STATISTICAL 

DESCRIPTION. 



The aim of this section is to connect the usual language of the ref. I^Sf model with the 
language of the statistical mechanics of the random media. 

Let us recall briefly the model of ( [|1^]). The model describes the QCD media as the 
usual Debye plasma, with a potential 

Vt{q) = gimcc expi-tqx,)/{f + i?) (Al) 

Here is the Debye screening mass. The parton-parton interactions are proportional to 
tr{T^Tj). The averaged potential vanishes everywhere since < V°' >^ TrT°' = 0, while the 
parton-parton interactions are proportional to 

Tr{T,^T!;) = 6ai,6,, {d^/dA)C2i (A2) 

Here di is the dimension of the SU(N) representation, C2i is the second Casimir, 
[[N'^—)/{2N) for quarks {di = N) and N for gluons {di = N"^ — 1). Thus after averag- 
ing we see that such model looks like multicomponent plasma, with a particular interacting 
charge for each set of color partons, with different types of partons not interacting in the 
leading order. Thus, we expect that apart from some numerical factors, connected with the 
Casimirs, the propagation of the wave packet in such a plasma will be qualitatively the same 
as it's propagation through the usual electromagnetic plasma described in the article. The 
true color effects will appear only when we shall take into account multiparticle correlations. 

We now turn to the description of the charged plasma in the language of the statistical 
mechanics ( pSf). In order to do it, we can just rewrite the statistical sum of the plasma as a 
Euclidean Functional Integral. The choice of the action must be such that it reproduces the 
correlation function ( |2.2| ). Clearly, the functional integral is Gaussian and can be represented 
as 

Z = JdE exp(- J d^qV{q^ + k^)V) (A3) 

Moving into the oscillator rest frame and taking into account that the time derivative term 
will be dominant, we come to the action ( |2.62| ). This description is just the usual language 
of the statistical mechanics of the random media. 
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APPENDIX B: LORENTZ TRANSFORMATIONS AND CENTER OF MASS 

MOTION. 



Let us introduce the center of mass coordinates (for simplicity let us assume that both 
particles have the same mass ) v and the relative motion coordinate u: 

v = {f^ + f2)/2 (Bl) 
M=(fi-f2)/2 (B2) 

(B3) 

We shall denote M the reduced mass of the oscillator in cm. reference frame. Then the 
Lagrangian can be written as a sum: 

L = j ds{v) + L,,,,{u) (B4) 

Indeed, the full Lagrangian is 

L = j ds{n) + ds{r2) - {V{n) - V{r2) + V{u))dt (B5) 

Here 



ds{r) = Jl-{^^ydt (B6) 



We then write: 



ds{n) = sl{l-{^^^±^Y)dt (B7) 



dsin) = p - {^-^^^Y)dt (B8) 

In these two expressions we expand their r.h.s. in Taylor series in u around v. Then we 
immediately obtain: 

ds{ri) + ds{r2) = ds{v) + [ dt{v)/{\Jl-u /2 = ds{v) + / dT{v)/2 (B9) 



Here T is the oscillator proper time. Adding the interaction term to the eq. ( |B9| ) we 
immediately obtain the decomposition (p^) . 

Let us now discuss briefly the Lorentz transformation [|17|. In the laboratory system we 



have only the electric field described by the static potential V. The Lorentz transformation 
says: 

V ={y' - vA'^)^ ^3 = {A'^ - vV')^ (BIO) 

In the laboratory frame A3 = 0, hence A'^ = vV . Substituting the latter formula in eq. |B1CI| 
we get 
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V = V'{1 - v^)-f = V'-f 



(Bll) 



Finally we get 



V' = 1/7 A' = vV 



(B12) 



Differentiating the latter equations in the moving reference frame we obtain the transforma- 
tion laws in the text. 



APPENDIX C: THE CALCULATION OF THE INTEGRAL. 

Here for completeness we shall write the formulae for the double integral used in the 



text: 



/infty POD 
dx / dyexpi-Qx'' - Q*y^ - 2Pxy)Hn{x)Hn{y) I {Tnl) 

( g + Q* + p2 - 1 - |g'|)"/V(IQ'l - p')("+^)/') 



Pn {-pi^{m - p'){Q +Q*- m +p'- 1)) (ci) 

The integral was calculated using ref. The Legendre Polynoms are defined as |l32| : 



2"n!P„(x) 



{x' - 1)" 



(C2) 
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FIGURES 
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FIG. 1. Probability of the oscillator that is a small wave packet with radius a=0.1 (here and 
below in the ground state radius units) to go from the ground to the tenth state as the function 
of time. We use here and below: = 1 GeV, M = 1 GeV-the reduced mass of the sample 
quarkonium, 7 = 30, and k = 0.12 GeV; = o.012 GeV^. 
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FIG. 2. Probability of the small oscillator wave packet with radius a=0.1 to be in the second 
state at the time t. 
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FIG. 3. Probability of the small oscillator wave packet with radius a=0.1 to be in the tenth 
state as the function of time. 
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FIG. 4. Probability of the oscillator to remain in the ground state, if it enters the media in the 
ground state, as a function of time. 
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FIG. 5. Probability of the oscillator to go from the ground to the second state as the function 
of time. 
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FIG. 6. Probability of the oscillator to go from the ground to the tenth state as the function 
of time. 
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FIG. 7. Probability of the charmonium entering the media in the ground state to remain 
charmonium at time t. 




200 400 600 800 1000 
FIG. 8. Density radius of ID oscillator in the random media as a function of time for big time 
scales, if initially it was in the ground state. Here and below Q{t) = x'^{t) 
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FIG. 9. Density radius of ID oscillator in the random media as a function of time for average 

time scales, if initially it was in the ground state. 



Q(t) 



0.25004 



0.25003 



0.25002 



0.25001 



.2 0.4 C.6 C.8 1 " 
FIG. 10. Density radius of ID oscillator in the random media as a function of time for short 
time scales, if initially it was in the ground state. 
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FIG. 12. Density radius of ID oscillator in the random media as a function of time for average 
time scales, if initially it was a small a=0.1 wave packet. 
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FIG. 13. Density radius of ID oscillator in the random media as a function of time for short 
time scales, if initially it was a small a=0.1 wave packet. 



31 



